Microscopic Theory for the Markovian Decay of Magnetization Fluctuations in 

Nanomagnets 



O 

o 

(N 

> 

o 
O 

ON 
(N 



Ioannis Rousochatzakia3 

Institut de theorie des phenomenes physiques, Ecole Polytechnique Federate de Lausanne, 
CTMC, BSP - Dorigny, CH-1015 Lausanne, Switzerland 

We present a microscopic theory for the phonon-driven decay of the magnetization fluctuations 
in a wide class of nanomagnets where the dominant energy is set by isotropic exchange and/or 
uniaxial anisotropy. Based on the Zwanzig-Mori projection formalism, the theory reveals that the 
magnetization fluctuations are governed by a single decay rate w c , which we further identify with the 
zero- frequency portion of the associated self-energy. This dynamical decoupling from the remaining 
slow degrees of freedom is attributed to a conservation law and the discreteness of the energy 
spectrum, and explains the omnipresent mono-exponential decay of the magnetization over several 
decades in time, as observed experimentally. A physically transparent analytical expression for u c is 
derived which highlights the three specific mechanisms of the slowing down effect which are known 
so far in nanomagnets. 

PACS numbers: 75.50.Xx, 75.75.+a, 76.60.Es 
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I. INTRODUCTION 

Over the past two decades there has been an intense 
theoretical and experimental interest in the field of molec- 
ular nanomagnets^ in view of their importance in the 
general context of quantum magnetism but also for po- 
tential applications, e.g. in quantum computing^, mem- 
ory storage^, and magnetic imaging 6 . The dominant en- 
ergy in these "zero-dimensional" magnets is set by in- 
tramolecular isotropic exchange interactions and/or easy- 
axis anisotropy. Two of the most central motifs in this re- 
search field have been£ the fact that the decay of the equi- 
librium fluctuations J(t) of the total magnetization S z is 
characterized by a single correlation frequency parame- 
ter lj c (or inverse decay time 1/t c ) over several decades 
in time, and the experimental manifestation of a dra- 
matic slowing down of lo c which establishes already from 
relatively high temperatures T. 

The correlation frequency parameter is a quantity of 
broad experimental relevance and of great theoretical 
value. For instance, according to linear response theory 
and the fluctuation-dissipation theorem, lu c is probed by 
ac (longitudinal) susceptibility and magnetization relax- 
ation measurements. It may also be measured directly 
by Electron Spin Resonance^ since uj c is the electronic 
spin-lattice relaxation rate (1/Ti) e i- Moreover, if uj c is in 
the regime of the nuclear Larmor frequency ojl, then it 
can be indirectly probed by Nuclear Magnetic Resonance 
(NMR) as in Antiferromagnetic Rings&2J£. Similarly, 
u>c is also of relevance in Muonii and Mossbauer— spec- 
troscopies. Previous theoretical studies of a;j2 i 13 i 14 i 15 i 16 
are based on the Master Equation (ME) which embodies 
the coupled dynamics of the full set of populations pi of 
the various eigenstates \l) of the isolated spin Hamilto- 
nian Ti. s . In this approach, one first evaluates the various 
phonon-driven quantum-mechanical transition rates Ww 
from \l) to |/') and then forms and diagonalizes numeri- 
cally the associated relaxation matrix 1Z. This procedure 



gives a multi-exponential decay. To account for the ob- 
served behavior it is then sometimes argued ^ 14 ' 15 ! 16 that 
uo c can be assigned to the lowest eigenvalue (or decay 
rate) Ao of 1Z. Despite being physically quite plausible, 
such a hypothesis alone does not guarantee that Ao cor- 
responds to the decay of the observable of interest (here 
S z ), since other observables (such as the magnetic energy 
?i 8 ) decay slowly as well. As shown by Santini et 
a full numerical treatment of this problem based on the 
ME approach requires also the calculation of the weight 
of each eigenvalue of 1Z in J(t). In agreement with ex- 
perimental data, the ME calculations of Santini et al~ 
showed that J(t) is dominated by a single eigenvalue and 
which, indeed, does not always equal Ao- A principal the- 
oretical question however still remains as to whether or 
not there is a physical reason which necessitates the ob- 
served dynamical decoupling of S z over several decades 
in time and in a large class of magnetic molecule clusters. 

In this article we put forward a theoretical framework 
which focuses right from the start on the principal quan- 
tity of interest, namely the fluctuations J(t) of the total 
magnetization S z . This approach reveals immediately 
the physical origin of the observed Markovian decay: The 
self-energy £(i) of the magnetization embodies fluctu- 
ations of non-conserved observables which decay much 
faster than S z . It is shown that this separation of time 
scales is guaranteed by the discreteness of the magnetic 
energy spectrum. The identification of the precise eigen- 
value of 1Z which governs the decay of J(t) with the zero- 
frequency portion of the self-energy and the derivation of 
an analytical expression for it (cf. Eq. I|27p or (pTlj) ') in 
lowest order in the magnetoelastic coupling is a second 
achievement of this approach. This expression is an op- 
timal starting point for direct and quantitative compar- 
ison with experimental data and thus allows to explore 
the specific low-energy phonon modes which are involved 
in the magnetoelastic process in all nanomagnets, an is- 
sue which is of broad interest and is currently largely 
unexplored. At the same time, this expression readily 
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demonstrates the three specific mechanisms of the slow- 
ing down effect which are known so far in nanomagnets: 
(i) the existence of an easy-axis anisotropy, (ii) the in- 
creasingly gapped character of the spectrum at low en- 
ergies, and (iii) the possibility of a crossover to a critical 
slowing down. 

The first is present in Single Molecule Magnets 
(SMM's)2 ' 17 i 18 which are perhaps the most intensively 
studied nanomagnets since, apart from the above "super- 
paramagnetic" behavior, they also provide a remarkable 
manifestation of magnetization tunneling^. Here, the 
presence of competing exchange interactions stabilizes a 
ferrimagnetic ground state of a large total spin S, e.g. 
S = 10 for the Mn 12 and Fe§ clusters. The energy scale 
within this large S manifold is set by an effective uniaxial 
anisotropy — \D\S^, giving rise to a correspondingly large 
barrier (U ~ 62 K for Mn 12 and 24.5 K for Fe 8 ^^ for 
thermal relaxation. The relaxation time t c of S z follows 
a phenomenological Arrhenius law r c oc e u ^ ksT which 
can reach values as long as 10 5 — f 7 sec at low T. This 
behavior is also revealed by complementary ac suscepti- 
bility measurements^. 

The second mechanism appears e.g. in Antiferromag- 
netic Rings AFMR's. These are planar ring-like magnetic 
clusters comprising an even number N of spins s whereby 
the dominant energy is set by the nearest-neighbor anti- 
ferromagnetic (AFM) exchange J. This stabilizes a non- 
magnetic S = ground state separated by a spin gap 
A ~ AJ/N (Lande interval rule) from the lowest S = 1 
excitation. Here there is no energy barrier and thus the 
corresponding uj c are on a much faster scale than the ones 
in SMM's. Still, it has been found£&±£ that oj c drops by 
several (three up to five) decades as T drops below A/ks- 
This dramatic decrease is systematically revealed by an 
enhancement of the nuclear spin-lattice relaxation rate 
1/Ti at T ~ A/fcs&SilA. A similar behavior of 1/Ti is 
observed in other AFM clusters but also in SMM's2£. In 
the latter, the analysis of 1/Ti is hindered by the NMR 
wipe-out effect^. 

The third mechanism of the slowing down effect 
emerges in Single Chain Magnets (SCM's)2ii22i22i2£i2£ 
which comprise finite ferromagnetic chains of variant 
lengths L. Here, a crossover from an activated to a criti- 
cal slowing down occurs as we tend toward the T = 
(thcrmodynamically) critical point, and in particular 
when the correlation length £ approaches L. Our ap- 
proach accommodates such a mechanism exactly as in 
critical bulk systems^ since it reveals that ui c is inversely 
proportional to the isothermal magnetic susceptibility xo 
(cf. Eq. (Zb). 

The remainder of this article is organized as follows. 
In Sec. [ii] we describe the main terms of the spin Hamil- 
tonian of the magnetic cluster and its interaction with 
the host lattice which drives the relaxational dynamics. 
In Sec. IIII Al we introduce the so-called Mori correla- 
tion function and express the main quantities of inter- 
est in terms of it. The Zwanzig-Mori projection tech- 
nique is then employed in Sec. IIII Bl in order to intro- 



duce the self-energy function and highlight its major role. 
A perturbative expansion of this quantity is then given 
up to lowest (second) order in the magnetoelastic cou- 
pling. In Sec. IIII CI we provide a rigorous justification of 
the Markov approximation based on the discreteness of 
the magnetic energy spectrum. We then arrive at our 
expression for the correlation frequency and provide a 
qualitative analysis of the three specific mechanisms of 
the slowing-down effect which are known so far in nano- 
magnets. In Sec. IIVI we revisit the problem of nuclear 
spin-lattice relaxation in AFMR's and clarify a number 
of aspects in the light of our treatment. Finally we give 
a short summary and overview of our results in Sec. |Vj 

II. HAMILTONIAN 

At not too low T, and given the typical energy scale of 
the magnetic excitations (a few tens of degrees Kelvin), 
the decay of the electronic spin fluctuations must origi- 
nate from the coupling to the deformations of the host 
lattice, and in particular, the low-energy phonon modes. 
The Hamiltonian of the full spin+phonon system can be 
generally written as 

H = n s + H ph + V , (1) 

where H s stands for the isolated spin Hamiltonian, H p h 
is the phonon energy and V denotes the relevant spin- 
phonon interaction channel(s). The first is taken to in- 
clude the most dominant magnetic energy terms 

% = J2 J v s i ■ S J " \ D \ S * + 9VbBS z , (2) 

i.e., the isotropic Heisenberg exchange interactions (lead- 
ing term), a possible effective easy-axis anisotropy (sec- 
ond term), and finally a Zeeman term for the coupling 
between an external longitudinal (along the easy z-axis) 
field B and the total magnetic moment S z = s l ■ As 
usually, g denotes the electronic spectroscopic factor and 
fiB is the Bohr magneton. In what follows, we shall de- 
note by \l) and E\ the eigenstates and eigenvalues of H s 
and by Mi the corresponding eigenvalues of S z . In our 
treatment we have excluded Dzyaloshinskii-Moriya inter- 
actions and transverse on-site anisotropics which are very 
small (< 1 K) and are relevant at ultra low temperatures 
only. As we show below, the observed mono-exponential 
decay of S z can be explained by the most dominant en- 
ergy terms which are included in Eq. (|2|). 

For the magnetoelastic coupling V, we may write with- 
out loss of generality 

V = J2 A 1 <8> B" , (3) 

where A 1 and are spin and phonon operators re- 
spectively. The latter represent low-energy strain fields 
probed by the electronic spin degrees of freedom. Details 
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about the o rigin a nd pos sible types of V are discussed 
e.g. in Refs. [I1MH28I2! As we show below, the major 



phonon quantities that govern the dissipative dynamics 
of the electronic spin fluctuations are the spectral densi- 
ties 



dt 



c (B^'(t)) 



(4) 



where the thermal average is taken over the canonical 
ensemble p ph (T) = e-^ H ^/Z ph , Z ph = Tr e ~ pn ^. 



III. MAGNETIZATION FLUCTUATIONS 

In the absence of spin-phonon coupling, S z is conserved 
since [S Z ,H S ] = 0. If the observables coupled to S z 
(through V) do not share this property it is physically 
expected that they decay faster than S z . This separation 
of time scales is at the origin of the Markovian decay of 
S z . And such a central physical ingredient can be ex- 
ploited naturally within the framework of Zwanzig-Mori 
projection formalis m 26 i 30 i 31 i 32 i 33 i 34 , whereby one may iso- 
late the important slow degrees of freedom by appropri- 
ately projecting out the remaining faster ones. To this 
end, it is expedient to first express the central quantities 
of interest in terms of the so-called Mori correlation func- 
tion which can be readily handled with the Zwanzig-Mori 
projection technique. 



Quantities of central interest in terms of the 
Mori correlation function 



Let us first introduce the Mori correlation function. 
We first note that any operator A can be represented 
as a "state vector" |A) in the Liouville space (space of 
operators). A convenient choice of a scalar product in 
this vector space is the so-called Mori product defined as 



1 r p 

(A\B) = - / dx(At Biihx)) = (A* 
P Jo 



2 _ e -pnc 



B), (5) 



where B(ihx) 



-hCx 



D 



-xH Be xH and £ = £ 



Cph + £v = £0 + £v denotes the Liouville operator. Its 
action on a given operator A is defined as 



in at 



= A 



(6) 



t=o 



Also, the thermal average is taken over the canonical en- 
semble p{T) of the full spin + phonon system 



p(T) = Z 



- y-K-m y - 



Z = Tr e 



-pn 



(7) 



where [3=1/ {kgT) denotes the inverse temperature and 
kg is the Bolzmann's constant. In what follows we as- 
sume that j3V <C 1 and thus the magnetoelastic coupling 



cannot influence static properties, i.e. the canonical den- 
sity matrix factorizes as 

p{T)~p s {T)®p ph {T) (8) 

where p p h(T) was defined previously while p s {T) — 
e~P n °/Z s , Z s =Tre _/3Ws . Now, the Mori correlation 
function M(t) is nothing else but the projection of the 
time-evolved "state" \5S z (t)) = e lCt \SS z ) onto \5S Z ), i.e. 



M(t) = (SS z \SS z (t)) = (SS, 



(3hC 



5S z (t)) , (9) 



where 5S Z = S Z - (S z ) and 5S z (t) = e mt / R SS z e- im ^ = 
e lCt 8S z . The one-sided Fourier (or Laplace) transform of 
M (t) reads 



M(s) = (SS z \SS z ) s = 



dt e tst (8S z \e lCt 8S z 



(10) 



On the other hand, the experimentally most relevant 
quantities which embody the dynamics of the magneti- 
zation fluctuations are the spectral density J{w) and the 
dynamical susceptibility x( s )- The first is the Fourier 
transform of the product correlation function J if), i.e. 



J{u>) 



dt 



l {5S z 8S z {t)) , 



(11) 



while the second is the one-sided Fourier transform of the 
commutator correlation function 



f 00 1 

X (s)= / dte lst (-[SS z ,SS z (t)}), s = u + ir), r)>0. 

JO 

The two quantities are connected via the fluctuation dis- 
sipation theorem 



2Ti 



(13) 



It is now straightforward to express the dynamical sus- 
ceptibility xi s ) m terms of A4(s) using their definitions 
Eqs. HUD and (TTJ]), and the Kubo identity (3(A\B) = 
— ( A[t4J, B] ) which follows from Eq. ((5]) by replacing B 
with B and using the relation (A t e~ /3n£ S) = (BA^). 
The resulting expression is 



X (s) = ~p(5S z \6S z ) s =l3(6S z 
= Xo + i(3s M(s) , 



£ 



£ 



SS Z 



(14) 



where xo is the isothermal differential susceptibility: 



Xo = P(6S Z \6S Z ) = p(6S; 
= p((6S z f) = - 



1 



D -(3hC 



1 d(S z ) 



-SS 2 



9Pb dB 



(15) 



The second line follows from /3£ ~ (3(C S + C p h) (since 
(3V < 1) and hL s S z = [H S ,S Z ] = 0. 
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B. Projection method and self-energy E(s) 

Having at hand Eqs. (fT4"| and (H2J), we may now con- 
centrate on the Mori correlation function which can be 
easily handled within the Zwanzig-Mori formalism. The 
idea is to project out from Eq. (|10|) fluctuations that cou- 
ple to S z through the resolvent Q = i(s + C)~ x . This can 
be done by employing the projection operators 



V 



\SS Z )(SS Z \ 

(ss z \ss z ) 



(16) 



and then by decomposing Q using the operator identity 
(C + V)- 1 = C- X -{C+V)- X VC~ X , withC = s + QC, and 
V = VC After some straightforward steps we obtain 



M(s) 



Xo 



P s + fl + zE(s) 



(17) 



where we have used M(t = 0) = (SS Z \8S Z ) = xo/P, and 
introduced the eigenfrequency f2 and the self-energy or 
memory function S(s) which are defined by 



iQ = — (dS z \i£5S z 
Xo 



S( S ) 



^(6S Z \CQ- 

Xo 



^(S Z \S Z ) 
Xo 

QC\SS Z 



QCQ 



(18) 
(19) 



The physical meaning of these quantities becomes trans- 
parent by taking the inverse Laplace transform of 
Eq. (fT7|) which reads 



We now employ a perturbation expansion for the eval- 
uation of the self-energy S(w + i0 + ). To this end, one 
needs to expand the resolvent Q = i(s + C)~ x of Eq. [21] 
in powers of Cy . Interestingly, S z is linear in V and thus 
we may already obtain the lowest order effect by simply 
replacing C by £q = C s + C p h in the denominator of 
Eq. (fT7|) . To second order then, (S Z \S Z ) w+i0+ equals to 



(S z 



C +uj + i0+ 



\S Z ) = (S z \tt6(jC +uj)\S z ) 



e f3hui _ i 

(S z S(£ +lj)S z ) , (23) 



/3hw 



where we have used the reality of ~E(uj + i0 + ). Employing 
the representation S(Cq + uj) = (2tt)^ 1 J dt e lCat , factor- 



iCnt 



izmg e 



Y,(uj+i0+) 



•iC s t AC-nht 



h , and using Eq. ([8]) we obtain 



P e 0Hw - 1 
Xo 2/3^ 



(24) 

where J v , (uj) was defined previously, while J", (uj) 
stand for the isolated spectral densities of the op- 
erators, i.e. 



(25) 



where the thermal average is taken over the canonical 
ensemble p s (T). 



iQM(t)- I dt"E(t')M(t-t'), i>0. (20) 
Jo 



This justifies the origin of the terms eigenfrequency and 
memory function for f2 and £ respectively. 

Equations (JT7J) , (JTSJ) and ([iU]) are exact. All informa- 
tion about the remaining fluctuations coupled to S z has 
been conveniently embedded into S(s). Here f2 — due 
to the definite time reversal signature of S z (also due to 
the Kubo identity given above). Importantly, this also 
implies that Qpz) = (1 — V)\S Z ) = \S Z ), which simpli- 
fies Eq. (PI) to 



E(«) - ^-(S z 
Xo 



+ 1~ Xo 



S z 



(21) 



Hence, E(cj + i0 + ) is proportional to the Mori correla- 
tion function (S^S^) tJ+i0+ and thus embodies the fluc- 
tuations of the "flux" operator— 



(22) 



where 
write th 
{l\^\V) 



rf- [S z , A* 1 ] . For later convenience let us 

iv — 



write the matrix elements of these operators: Tf,, 



5M n Af„, where 5M», 



Mil 



M,. 



C. Markovian behavior and correlation frequency 

The origin of the mono-exponential decay of S z follows 
immediately from Eq. (|24|) and this is related to the dis- 
crete character of the spectrum. The spectral densities 
J s (uj') of F/j, are peaked at the characteristic Bohr fre- 
quencies uj' = (Ei — Eii)/h of 7i s . On the other hand, we 
are interested in a long time regime: For all relevant ex- 
perimental frequency scales (e.g. ujl in NMR) uj ~ uj c <C 
uj' , and thus, given the extremely short phonon correla- 
tion times, we may replace^ J ph (uj — uj') ~ J ph (—uj'), 
which is equivalent with the usual Markov approxima- 
tion of replacing S(w + iQ+) ~ S(iQ+) in Eq. dTTJ) . This 
"memory-less" character of the self-energy means that 
the fluctuations of the observables coupled to S z (i.e. 
the observables J-^ that enter in S z ) decay in times t' c 
much shorter than the decay time r c = 1/uj c of M(t). 
The physical origin of this separation of time scales can 
be attributed to the fact that, in contrast to J 7 *, S z 
is a conserved variable for the isolated magnetic clus- 
ter and thus is expected to decay much slower in the 
presence of the small magnetoelastic coupling. At the 
same time, this separation of time scales is at the ori- 
gin of the Markovian decay of S z : As long as t 3> t' c , 
one may replace M(t — t') with M(t) inside the integral 
of Eq. (|20[) and extend the upper limit to infinity, i.e. 
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f- t M(t) ~ - J °° dt'E(t') ■ M(t) = -S(iO+) • M(t), which 



give; 



,3!) 



M(t) = (xo/P) e~^ + )\ t > 



(26) 



We can thus identify the zero-frequency portion of the 
self-energy with the correlation frequency lu c : 



^AO , J —oo 



On the other hand, according to Eqs. (fl4|) and (|13j) . the 
susceptibility and the spectral density have the familiar 
Debye (or Lorentzian) form: 



JM = 



S + Vj) c 
p—f3hoj 



(28) 
(29) 



Equations (|2~7| , (|28p and (f2"9"|) are our central results. It is 
expedient to rewrite Eq. (|27p in a more convenient form, 
in terms of the various quantum-mechanical transition 
rates Ww from \l) to \l'). To this end, we use the spectral 

representation J^,{uj') = 2ttJ2w Pi^w^vi^' ~ 
where Pi = e~^ El /Z s and !F W were given previously We 
find 



2tt o.^ 1 



LI' 



MM' 



(30) 



The last factor of this equation can be identified to be 
the quantum-mechanical transition rate Ww from \l) to 
\l'). We therefore obtain the convenient expression 



2tt 



1,1' 



s -^w u , 

(SSI) 



(31) 



We remark here that we can eliminate the factor of 1/2 
appearing in the right hand sides of Eqs. (127]) and (|3"Tj) by 
restricting the integral of (f27|) to positive (or negative) 
frequencies only, or by constraining the sum of pip to 
run over I < I'. This is allowed by the detailed balance 
condition J^,{-oj') = e _ ^ W J^(uj') (i = s, ph) or its 
equivalent form PiWw — PyWin. 

Let us now describe in qualitative terms the main fea- 
tures of Eq. (f27|) or Eq. ([31]) . According to these expres- 
sions, uj c is a weighted sum over all relevant resonance 
channels between the spin and the phonon system. Each 
contribution enters in a transparent and physically ap- 
pealing factorized form. Energy conservation is ensured 
by the opposite frequency arguments in J s and J ph . 
The one-phonon portion of the latter is proportional to 
the number of available (i.e. thermally excited) phonons 
in resonance with the electronic spin system. We also 
note that it is and £> M triggering the decay of S z (cf. 
(gl]) ) that enter in Eq. (gTJ). This is also reflected in 
Eq. PT]) by the fact that transitions between two levels 



1,1' contribute to uj c only if they differ in their magnetic 
moment. More generally, each of the three known mecha- 
nisms of the slowing down effect can be now associated to 
one of the three major quantities appearing in Eq. (|27p . 
namely J ph , J s and xo- In AFMR's (see also Sec. [TV] 
below), it is the drastic reduction of the resonant por- 
tions of the phonon spectral weight J ph as T drops be- 
low the spin gap. In SMM's on the other hand, the effect 
appears already well above the typical intra-barrier ex- 
citations, and thus must be attributed to selection rules 
in J s which necessitate an "over the barrier" relaxation 
process. Finally, the appearance of the susceptibility in 
the denominator of Eq. (|27p is underlying the critical 
slowing down observed in SCM' o 21 i 22 i 23 i 24 i 25 similarly to 
bulk ferromagnets^. 

We should finally note here that one may start with 
the magnetic energy 7i s instead of S z and show, follow- 
ing exactly the same steps as above, that this observable 
decays also independently from the remaining slow de- 
grees of freedom. Its corresponding decay rate u>e can be 
also derived and the result is fully analogous to Eq. (|3ip : 



UJE_ 

2^ 



1,1' 



l (5Hl) 



(32) 



where 5E vi = E v - E h {5U 2 S ) = k B T 2 C m , and C m is 
the magnetic specific heat. The "magnetic energy cor- 
relation frequency" loe is of experimental relevance in 
time-dependent specific heat measurements^ and is dis- 
tinct from lo c . It is nevertheless expected, based on the 
above discussion, that oje also manifests a dramatic slow- 
ing down similar to lu c . 



IV. NUCLEAR SPIN-LATTICE RELAXATION 
RATE IN AFMR'S 

Here we revisit the problem of proton 1 jT\ in AFMR's 
and clarify a number of aspects in the light of Eqs. ([29]) 
and (|2"T|) or ([3T|) . Quite generally, since the proton spin 
is a local probe, 1/Ti samples all (auto- and pair-) elec- 
tronic spin spectral densities evaluated at u>l. Remark- 
ably however, in clusters with equivalent ionic spins (such 
as AFMR's) symmetry arguments together with the fact 
that lol is much smaller than the typical Bohr frequen- 
cies u>b assert that away from level crossings the relevant 
terms contributing to 1/Ti are all proportional to the 
spectral density of the total magnetization^^, i.e. 



= AJ(u L ) , 



(33) 



where the proportionality constant A embodies the spe- 
cific details regarding the geometry and the strength of 
the nuclear-electron hyperfine coupling- . The remaining 
terms of the electronic spin spectral densities are van- 
ishingly small and become relevant only at field-induced 
level crossings^. For symmetric clusters then and using 
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Eq. {29]) (with [3huj L <C 1) we obtain 

k = ^ a ^ t ^tBty • ^ 

which is the Lorentzian form used previously 8,9 . 

In AFMR's, J s (u b )/(xqT) does not change dras- 
tically over the whole temperature range of interest. 
Hence, the observed dramatic decrease of u> c and the 
resulting enhancement 8,9,10 of 1/Ti when w c (T) = u) L 
(cf. Eq. (|34p) stems, as mentioned above, from a dras- 
tic reduction of the number of available (thermally ex- 
cited) phonons which are in resonance with the magnetic 
cluster. This is because on decreasing temperature the 
relevant low-lying magnetic excitations are increasingly 
gapped while the frequency peak of J ph {ui) shifts toward 
lower uj. Eventually the phonon spectral density peak 
shifts below A/h and one-phonon processes become in- 
effective in driving the electronic spin fluctuations, thus 
resulting in a drastic drop of u> c . One can also explain the 
observed field independence of w^iiS at intermediate 
temperatures and for B <C A/ (5/15) as follows. The rel- 
evant Bohr frequencies (with 5Mpi ^ 0) in these temper- 
ature and field regimes are A/?i±a; e uj e and 2u> e (here u e 
is the electron Larmor frequency) . These frequencies cor- 
respond respectively to the transitions between the lowest 
singlet and the first triplet excitation (with SMyi = ±7i) 
and to transitions within the lowest triplet. The later do 
not contribute appreciably since J ph (Lo e ) <C J ph (A/K). 
At the same time J ph (A/h ± u e ) ~ J ph (A/h) and 
Xo{B) — Xo(0). As a result uj c is field-independent at 
intermediate temperatures and B <C A/^g/is)- A field 
dependence should only arise at higher fields and very 
low temperatures or for clusters with small A. 

V. SUMMARY 

We have presented a microscopic theory for the om- 
nipresent mono-exponential slowing down effect observed 



in a wide range of molecular nanomagnets. To this end, 
we first expressed the magnetization fluctuations in terms 
of the associated self-energy or memory function which 
embodies the fluctuations of the observables that couple 
to S z . It is then shown that the discreteness of the mag- 
netic energy spectrum guarantees that the self-energy de- 
cays in times much shorter than l/u> c . This separation of 
time scales results in a dynamical decoupling of the mag- 
netization fluctuations and thus explains their observed 
Markovian behavior. In addition, we have derived an an- 
alytical expression for the correlation frequency parame- 
ter lo c which embodies in a convenient factorized form all 
essential ingredients. This expression highlights the three 
specific mechanisms of the slowing down effect which are 
known so far in nanomagnets, namely (i) the existence of 
an anisotropy barrier, (ii) the increasingly gapped char- 
acter of the magnetic energy spectrum at low energies, 
and (iii) the possibility of a critical slowing down. In ad- 
dition, this formula contains, in a convenient factorized 
form, the noise spectra of the relevant low-energy strain 
fields which are at present largely unexplored. Hence this 
work paves the way for a more direct and quantitative 
comparison with experiment and a deeper understanding 
of the underlying microscopic relaxation channels present 
in all nanomagnetic systems. 
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